Global defects in field theory with applications to condensed matter 
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We review investigations on defects in systems described by real scalar fields in (D, 1) space-time 
dimensions. We first work in one spatial dimension, with models described by one and two real 
scalar fields, and in higher dimensions. We show that when the potential assumes specific form, 
there are models which support stable global defects for D arbitrary. We also show how to find 
first-order differential equations that solve the equations of motion, and how to solve models in D 
dimensions via soluble problems in D = 1. We illustrate the procedure examining specific models 
and showing how they may be used in applications in different contexts in condensed matter physics, 
£NJ ■ and in other areas. 



■ I. INTRODUCTION 

The search for defect structures of topological nature is of direct interest to high energy physics[l-3] and condensed 
matter [4, 5]. In the present work, we shall be mainly interested in defect structures, as they appear in high energy 
physics. However, we shall investigate models which may be of good use in applications in condensed matter, mainly 
in pattern formation and in magnetic systems, where one finds rich scenarios with defect structures playing important 
role to control some of their physical features. 

In Field Theory, the search for defect structures started several years ago. In the 1960's, works of Hobard[6] and 
Derrick[7] have shown that for standard real scalar field models, there is no way of generating stable defect structures 
£j ' unless one works in (1, 1) space-time dimensions. In one space dimension, the defect structure of topological nature 
is named kink. This result has led to the standard route to defect solutions, in which one enriches the model 
to support vortices, if one works with (2, 1) space-time dimensions, with the introduction of Abelian gauge field, 
coupled to charged scalar matter[8]. Or to support monopoles, if one works in (3, 1) space-time dimensions, with the 
[ — . introduction of non- Abelian gauge fields coupled to charged scalar matter [9-1 1]. 

' For models described by a single real scalar field, the Lagrange density has the form 
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We work in (D, 1) space-time dimensions, and the metric is (+, — ), with x^ = (xq = x = t, x\ = —x , X2 

—x 2 , ...,Xr> = —x D ). The equation of motion is 



d 2 S , dV 



For static field configuration, the equation of motion becomes 



(3) 



The energy can be written as 

dx 1 dx 2 ...dx D 



(4) 



In this work we shall be interested in models described by scalar fields, defined in D = 1 and in D > 2. The case of 
a single spatial dimension is standard and will be investigated in the next Sect. II. Extensions to two or more spatial 
dimensions will be considered in Sect. III. 

In D = I we consider models described by one field 

v(4>) = \ K (5) 
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and by two fields[12-14] 

V(<f>,x) = \(W$ + W*) (6) 

Here W is a smooth function of the fields, and W$ — dW/d<p, etc. For higher spatial dimensions, for D > 2 we work 
with models described by the potentials[15] 

n<i>;x) = \^wl (7) 

and 

V(4>,x;x) = l^ (wl + wl) (8) 

The present generalization to higher dimensions is different from the extensions one usually considers to evade the 
Derrick-Hobard theorem, which include for instance constraints in the scalar fields and/or the presence of fields with 
nonzero spin - see, e.g., Ref. [3], and other specific works on the subject[16, 17]. Potentials of the above form appear 
for instance in the Gross-Pitaevski equation, which finds applications in several branches of physics[18]. 



II. ONE SPATIAL DIMENSION 



Here we work in (1, 1) space-time dimensions. We investigate models defined by potentials that can be written as 
in Eqs. (5)-(6). Our main motivation is to show specific features of models described by one and by two real scalar 
fields. 



A. One Field 



For models described by a single field, the equation of motion for static solutions has the form 



We are searching for solutions that obey the boundary conditions 



d<f> 



lim 4>(x) — ► (f>, lim ► 



(9) 



(10) 



where the singular point <j> of the potential obeys V{4>) = 0. 

For solutions which obey the above boundary conditions, the second-order equation of motion is equivalent to the 
two first-order equations [19, 20] 



d0 
dx 



= ±W 6 



(11) 



These first-order equations also appear in the investigation of the energy of static solutions. The energy density has 
the form 



s{x) 



1 



dx 



(12) 



It can be written as 



. , dW 1 fd<f> „, V 



(13) 



The energy is obtained after integrating the energy density in the full real line. It is positive, and it can be minimized 
to the value E — \ AW\ — \W[<j>(oo)] — W{4>(— oo)]| for field configurations that solve the first-order equations (11). 
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Solutions of the first-order equations (11) are named BPS and antiBPS states, and are classically or linearly stable. 
We show this by considering <j>{x, t) = cj>(x)+r)(x, t) and using the equation of motion to get, after discarding non-linear 
contributions in 77, 



d 2 r] d 2 rj d 2 V 



+ 



r] = (14) 



dt 2 dx 2 d(f> 2 
The static field <p(x) is time- independent, then we can use 

v( x , t) = ^2 ^(x) cos ( w «*) ( 15 ) 

n 

to obtain 

d 2 rj n 



dx 2 



+ U(x)rj n (x) = w n r) n (x) (16) 



where U(x) = d 2 V/d<fi 2 has to be calculated at the classical static configuration <t>{x) that describes the defect solution. 
The equation (16) is a Schrodinger-like equation, which leads us to Quantum Mechanics. The Hamiltonian is 

H = -^ + U(x) (17) 
For the potential V = (1/2) Wj, it can be written as 

H=-^ + W 2 , + W^W 4>4>4> (18) 
The structure of the field-theoretical model furnishes a pair of Hamiltonians, H± = S±S±, defined by 

S ± = ±±W„ (19) 

which nicely lead to supersymmetric quantum mechanics [21, 22]. For the issue of stability, the important point is 
that both H± are non-negative, and so no eigenvalue in Eq. (16) can be negative, ensuring linear stability of the BPS 
states. 

To illustrate the general situation, we consider the model defined by the potential 

V{4>) = \\\v 2 -<f) (20) 
where A and v are parameters, real and positive. This is the (j) 4 model, and here W has the form 

W((f>) = \v 2 <j>- \\<j> 3 (21) 

with Wcf, = \(v 2 — (j) 2 ). The potential has two minima, at <f>± = ±v. The energy of the static solutions has the form 
E = (4/3)Aw 3 . 

The model admits static solutions in the form 



4>± {x) = ±v tanh(Aux) (22) 

where we are choosing x — as the center of the defect. These solutions are named kink (+) and anti-kink (— ). They 
have two main features: the amplitude (v) and the width (1/Xv). Kinks and anti-kinks can be embedded in (2, 1) 
space-time dimensions, to give rise to domain ribbons, or in (3, 1) dimensions, to generate domain walls. In the above 
example, the defect is structureless and corresponds to the Ising wall one finds in magnetic systems [4]. 

In the (f> 4 model, the two minima defines a single topological sector. Both minima are asymmetric, that is, they 
break the Z 2 symmetry of the model. The choice of the minimum leads to the asymmetric or ordered phase; thus, if 
we introduce a mechanism to restore the symmetry, we should necessarily describe a second-order phase transition. 

Another example, similar to the above one is given by the potential 



V(4) = ±\ 2 (v-W) 2 



(23) 
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where v and A are positive. We name this the modified <p 2 model. It was used in[23] as an exactly integrable model, 
to illustrate the presence of a single-soliton solution in a system of discrete chains. It can be written in terms of 



The first-order equations are 



The BPS solutions are 



W = \v<j) - -\\<j>\<i) 



ax 



cj)(x) = ±2v 



tanh(^) 
1 + tanh ( ^ 



(24) 



(25) 



(26) 



There is a single topological sector, which supports BPS solutions. The energy of the BPS sector is given by E = Xv 2 . 
In Fig. [1] and [2] we plot potentials and defect solutions for both the </> 4 and modified 4> 2 models. 




FIG. 1: Plots of the potentials for the 4> 4 an d modified 
v — 1 and A = 1, for simplicity. 



models, depicted with solid and dashed lines, respectively. We use 



Another example is given by the potential 



V(0) = -AV(u 



A2^2 



(27) 



where v and A are positive. This is the <fi 6 model, and it has three minima, one at tfto = 0, which preserves the Z% 
symmetry of the model, and two at <fi± = ±v, which break the Z2 symmetry. Here we have W — Av 2 <fr 2 /2 — A^> 4 /4. 
The energy of the BPS solutions is E = (l/4)Ai> 4 . The solutions are 



'j|(x) = ±«-\/ ^[1 ± tanh(Awa;)] 



(28) 



The <p 6 model supports two topological sectors, degenerate in energy, defined by the minima —1 and 0, and and 1. 
Furthermore, the minimum at zero does not break the symmetry of the model, so it represents a symmetric phase. 
Differently from the <p model, the </> 6 model supports two distinct but degenerate phases, one symmetric and the 
other asymmetric; thus, it can be used to describe a first-order phase transition. 

For this reason, it seems to be interesting to describe another model, similar to the above one, but with lower power 
in the field. Such a model was invented very recently in [24], and it is governed by the potential 



1 



(29) 



FIG. 2: Plots of the kink-like solutions for the <fi 4 and modified 2 models, depicted with solid and dashed lines, respectively. 
We use v = 1 and A = 1, for simplicity. 



where v and A are positive. We name this the modified <fi 4 model, and it can be written in terms of W = \v(f> 2 /2 — 
X\(j)\(j) 2 /3. It also has three minima, one at </> = and two at (f> — ±v. However, the topological solutions are different, 
given by 



They have energy given by E = (1/6) Aw 3 . 

In Fig. [3] and [4] we plot the potentials and defect solutions for both the 4> 6 and modified (f> 4 models, to show how 
they differ from each other. There we see that the modified <fi model is more symmetric than the former <fi 6 model, 
and may be used in applications in condensed matter as symmetrical interfaces or fronts, which can also be of interest 
to model periodically forced oscilatory systems governed by the Ginzburg-Landau equation[25]. 



FIG. 3: Plots of the potentials for the (f> 6 and modified cf) 4 models, depicted with solid and dashed lines, respectively. We use 
v — 1 and A = 1, for simplicity. 

Another model is the double-sine- Gordon model. It is defined by[26] 




(30) 
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}v(*) 




-1 



V{4>) 
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[r + cos(</>)] 2 



(31) 
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FIG. 4: Plots of the kinks for the (f> 6 and modified </> 4 models, depicted with solid and dashed lines, respectively. We use v = 1 
and A = 1, for simplicity. 



where r is a real parameter; we consider r S (0, 1), and the limits r — > and r -> 1 lead to distinct sine-Gordon 
models. This potential can be written in terms of 



W{cj>) = 



vT+7' 



\rm + sm 



(<f>)} 



The first-order equations are given by 



[r + cos(0)] 



There are two kind of solutions, representing the large and small topological sectors; there are large kinks 



4>i(x) — ±2arctan 



(y^t,anh(VT-^)) 



and small kinks 



Their energies are, respectively, 



and 



4>s(x) = ±7r — 2 arctan 



yiZZtanh 7a;) 



Ei = 4\/l - r + 4r 



7r — arccos(r) 
arccos(r) 



^ = 4\/l - r - 4r 

vT+7 

Another example is given by the family of models introduced in [15] 

1 / , p-i , p+i \ 2 



(32) 



(33) 



(34) 



(35) 



(36) 



(37) 



(38) 



Here both the field and coordinates are dimcnsionlcss. The parameter p is real, and it is related to the way the field 
self-interacts. These models are well-defined for p odd, p — 1,3, 5... For p = 1 we get to the standard <p A theory. For 
p = 3, 5, ... we have new models, presenting potentials which support minima at <p = and ±1. For p odd the classical 
bosonic masses at the asymmetric minima <f> = ±1 are given by m 2 = 4/p 2 . For p — 3, 5, ... another minimum appear 
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at 4> — 0- However, the classical mass at this symmetric minimum diverges, signaling that = does not define a 
true perturbative ground state for the system. 
We consider p odd. The first-order equations are 



which have solutions 



d<f> ,£=i ,e±1 
—— = ±(p p ^ (p p 
ax 



4> { l' p \x) = ±tanh p (z/p) 



(39) 



(40) 



We consider the center of the defect at x = 0, for simplicity. Their energies are given by e[ p I = Ap/(4p 2 — 1), and we 
plot some of them in Fig. [5] . 



FIG. 5: Plots of kinks for p — 1, 3, 5, depicted with solid, dashed, and dot-dashed lines, respectively. 

We see that solutions for p — 3,5, connect the minima <j> = ±1, passing through the symmetric minimum at cj> = 
with vanishing derivative. They are new structures, which solve first-order equations, and we call them 2-kink defects 
since they seem to be composed of two standard kinks, symmetrically separated by a distance which is proportional 
to p, the parameter that specifies the potential. To see this, we notice that the zero modes are given by 

= c (1 ' p) tanf-^x/p) sech 2 (x/p) (41) 

where c {1 ^ = [{4p 2 - l)/4p] 1 /2 

- we use rig D ' p ' > to represent the zero modes. These zero modes concentrate around 
two symmetric points, which identify each one of the two standard kinks. The defect solutions are non-perturbative, 
and for p = 3, 5, ... they somehow know that the minimum at zero cannot be used as a perturbative ground state. 
The cases with p even are different. The case p — 2 is special: it gives the potential 

V{^)= 1 -<j ) -^+ 1 -f (42) 

which supports the non topological or lumplike solution 

; (1 ' 2) (x) = tanh 2 (x/2) (43) 

The lumplike solution is unstable, as we can see from the zero mode, which is proportional to tanh(a;/2)sech 2 (cc/2): 
the zero mode has a node, so there must be a lower (negative energy) eigenvalue. In the present case, the tachyonic 
eigenvalue can be calculated exactly since the quantum-mechanical potential associated with stability of the lumplike 
configuration has the form 



U(x) = l - 3sech 2 (x/2) 



(44) 
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This potential supports three bound states, the first being a tachyonic eigenfunction with eigenvalue Wq = —5/4, the 
second the zero mode, w\ — 0, and the third a positive energy bound state with w\ =3/4. 

The other cases for p even are p = 4,6, ... These cases require that <j> > in Eq. (38), but we can also change <f> — > — cj) 
in Eq. (38) and consider 4> < 0. We investigate the case <j> £ [0,oo); reflection symmetry leads to the other case. We 
notice that the origin is also a minimum, with null derivative. These models also support topological defects, in the 
form (f>( 1,p \x) — tanh p (x/p) (x > 0), with energies e[ p ], = 2p/(Ap 2 — 1) for p = 4, 6, ... These solutions solve the first 
order equation d<fi/dx = W^; the other equation d<p/dx — —W^ is solved by 4>(x) = — tanh p (x / p) (x < 0). 

In (1, 1) space-time dimensions, the case of p — 3, 5, ... are very interesting since they give rise to 2-kink solutions. 
The energy density of the 2-kink solutions are given by 

s(x) = temh 2p - 2 {x/p) sech 4 (x/p) (45) 

We plot this energy density in Fig. [6]; it is zero at the center of the defect, and it has maxima at the two symmetric 
points x± = ±p arcsechy / 2/(p + 1), showing that the 2-kink solution engenders internal structure, being composed 
by two kink solutions. It is very interesting to notice that this kind of behavior was recently found in a magnetic 
system, if one constrains the geometry of the material in a very specific way [27]. 

t s(x) 

0.1 




FIG. 6: Plots of the energy density for p — 1,3,5, with solid, dashed, and dot-dashed lines, respectively. The case p = 1 is 
depicted with the factor 1/10, to easy comparison. 



B. Two Fields 



For models described by two fields, the potential is given as in Eq. (6). The equations of motion for static solutions 
have the form 

= W+WM + WxWrt (46) 
= W^ X + W X W XX (47) 
The energy density for static solutions can be written as 

The energy is then minimized to the value E = |AW|, where 

AW = W[0(oo),x(oo)] - W[0(-oo),x(-oo)] (49) 
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for field configurations which obey the first-order equations 



We notice that the first-order equations can be seen as dynamical systems. In this case, the singular points solve 
W,p = and W x = 0, and so they are the absolute minima of the potential. We identify the topological sectors of the 
model by choosing pair of absolute minima: the pair (4>i,Xi) an d (4>j,Xj) identify the sector (ij), with energy 
If W(4>i,Xi) ^ W((j)j,Xj), the sector (ij) is a BPS sector, with energy given by — \W(4>i,Xi) — W(4>j,Xj)\- In 

this case, the first-order equations (50) have non-trivial solutions which connect the minima (4>i,Xi) an d {<f>j>Xj)- 

To exemplify the general situation, let us consider the model defined by 

W = \v 2 <P - ±A0 3 - M X 2 (51) 

This model was first investigated in[12], and it can be seen as an extension of the <j) 4 model considered in the former 
subsection. The first-order equations are 

g = ±A(t, 2 -^) T MX 2 , fr = T2/ ^ x (52) 

For A//i > 0, the absolute minima are v\_2 — (±v,0) and v^,a — (0, ±V\f\fjl). There are five BPS sectors, identified 
by the energies: E^ = (4/3)|Av 3 | and E (13 1 = E^ = E (23 1 = E (24 1 = (2/3)|Au 3 |. Also, there is a single non-BPS 
sector, identified by the minima and V4. 

In the sector defined by the minima v\ t 2 we have several BPS solutions. Some of them are given by 

4>±{ x ) = ±wtanh(Avx), x( x ) = (53) 

They describe a straight line segment connecting the two minima (±v, 0) in the ((f), x) plane. Other pairs of solutions 
are 

4>±{ x ) — ±wtanh(2/it>a;), x( x ) — ±«\/ A//x — 2 sech(2fivx) (54) 

for X/fi > 2. They describe eliptycal arcs connecting the minima V1.2 — (±v,0). These solutions describe orbits in the 
configuration space which obey 

' X 2 =v 2 (55) 



X/H - 2 



We notice that in configuration space, the one field solutions describe straight line orbits, and the two field solutions 
describe eliptycal orbits. They very much remind us of linearly and eliptycally polarized light. This scenario is 
interesting, and may be used to model more sophisticated problems. A good example can be found in works used to 
describe chiral interface[40, 41]. 

The presence of the second field x leads to the interesting situation: For x = 0, in the center of the defect the second 
field does not vanishes, and may be used to describe internal structure. The profile of the wall the model generates 
is the same of finds in the so called Bloch wall scenario, which appear in models described by the Ginzburg-Landau 
equation, as one finds in magnetic systems [4] and in other situations [5]. 

Models described by two or more fields are also used to describe conformational structures in almost one dimensional 
polimeric chains. Work in this line has been done in several different contexts, in particular to describe defect structure 
in a ferroelectric crystal[42, 43] and in polyethylene [44, 45]. Other form of interactions can be of use to describe 
conductivity in Langmuir films [46]. 

The presence of internal structure can be extended to describe scenarios where a single field entraps two fields, 
generating more sophisticated structures, where networks of domain walls live inside a given domain wall. This idea 
was first used to simulate nanotubes, fullerenes and soliton stars[47-49]. 



III. TWO OR MORE SPATIAL DIMENSIONS 



To find defect structures in higher dimensions, we consider the class of potentials given by Eq. (7). In this case, we 
suppose the model supports static and radial solutions, in the form <f> — 4>(r). This leads to the equation of motion 



d ( D _ 1 d<j)\ 1 



r 



dr V dr 



W+Wm (56) 
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which can be written as 

-l d ( r>-i# 



^ ^=]-^ (57) 



We change r — ► y, obeying r D 1 (d/dr) = (d/dy) to get 

dy 2 



= W*W^ (58) 



which is the equation of motion for static solutions in D = 1. Thus, we reduce the problem of finding static solutions of 
some (D, 1) dimensional model, to that of effectively finding static solutions of some much simpler (1, 1) dimensional 
model. 

The energy of the static, spherically symmetric field configuration has the form 



dr 



(59) 



where flu is the angular factor: Q D = 2tt d / 2 T(D/2), with T standing for the gamma function. We can work with the 
energy to write 



It can be minimized to the value E BPS = Q, D \AW\, where AW — W((f>(r — » oo)) — W((j)(r — > 0)), for field configura- 
tions that obey 

The above results are extensions to D > 2 of results valid in (1, 1) space-time dimensions. To get to them, an 
important manipulation was the change of variables r — > y; for D = 2 it leads to y = ±ln(r), which relates the 
interval [0, oo) to (—00,00). However, for D > 3, one gets y = ±r 2 ~ D /(2 — D), and now y must be in the interval 
[0,oo) or (—00,0]. Thus, for D > 3, the unidimensional problem that one has to solve is different from the standard 
situation, since the coordinate y does not span the full real line anymore. For this reason, to work in higher (D > 3) 
dimensions, one has to invent new models, which need to be solved in the interval [0, 00). 

The presence of first-order equations corresponds to the Bogomol'nyi bound, now generalized to models for scalar 
field that live in D spatial dimensions. The result shows that the solutions obey first order equations and have energy 
evenly distributed into gradient and potential portions. We notice that for D = 2 the model engenders scale invariance. 
Evidently, the solutions of the above first order equations solve the equation of motion (57) for potentials given by 
Eq. (7). Also, we follow former work[19] and introduce the ratio R = r D ~ 1 (d<fr/dr)/W ( f > . For field configuration that 
obey 0(0) = <j> and lim r ^o defi/dr one can show that solutions of the equation of motion also solve the first-order 
equations (61). This extends former result[19] to the present investigation: it shows that the equation of motion (57) 
completely factorizes into the two first order equations (61). Thus, all the topological solutions are BPS states. 

The equation of motion for <f> = 4>{r) is 



-W^W^ (62) 



dr 

The first order equation is given by (61), and their solutions solve the equation of motion and are stable against 
radial, time-dependent fluctuations. To see this explicitly, we consider (j>(r,t) = 4>(r) +J2k 1 lk(r)cos(wkt). For small 
fluctuations we get Hrjk — w\r\k, where the Hamiltonian can be written as 

H = ^=5 {-^l T w **) { rD "i T w **) (63) 

It is non-negative, and the lowest bound state is the zero mode, which obeys r D ~ x dr\{)l dr = ±W^r?o- This gives 
i]o(r) = c cxp (± J drr x ~ D W^) , where c is the normalization constant, which usually exists only for one of the two 
sign possibilities. We can also write 770 (r) = cM 7 ^, which is another way to write the zero mode. 



11 



In D = 1 one usually introduces the conserved current j M = e^d^cf). We see that p = dcj)/dx and so p 2 gives the 
energy density of the field configuration [31]. Thus, we introduce 



Qi 



dxp 2 



(64) 



as the topological charge, which is exactly the total energy of the solution. 

We now consider D = 2. The equation of motion is V 2 (/> = (l/r 2 )W c f,W^ ( p. We search for solution 4>(r), which only 
depends on the radial coordinate, obeying (f)(0) = 4> and lim r _,o d<p/dr — ^ 0. In this case we get 



dr \ dr 



dW d 2 W 



(65) 



We write dx = ±r~ 1 dr to get d 2 <f)/dx 2 = W^W^. This result maps the D = 2 model into the D = 1 model. 

Wc sec that r = cxp(±x), which shows that the full line x £ (-co, oo) is mapped to r £ [0, oo). We notice that since 
the center of the defect is arbitrary, so it is the point r = 1 in D = 2; thus the solution introduces no fundamental 
scale, in accordance with the scale symmetry that the model engenders. If one uses the model (38) with p odd to 
define the potential in this case, we get the solutions 



4>%*\r)=± 



r z/p _ i 
t 2 /p + 1 



(66) 



Their energies are = 8irp/(4p 2 — 1). In Fig. [7] we depict the defect solution for p = 1. The corresponding zero 
mode is given by rj^ '^(r) = y / 32~77r[r 2 / '(r 2 + l) 2 ]. This zero mode binds around the circle where the defect solution 



vanishes, as we show in Fig. [7], where we also plot p^' X \r) = 32r i /ir(r 2 + 1) 



-0.5 




FIG. 7: Radial defect solution cj)(r) and density p(r) of the corresponding zero mode for D = 2, in the case p = 1. 



We now investigate a model described by a complex scalar field. Here the equation of motion has the form 

1 dV 

r* dip 



(67) 



where bar indicates complex conjugation. We suppose that the field configuration has the form ip(r,9) — 4>(r)e lra , 
for m — 1, 2, 3, ... In this case we get 



r d_ (d£\ dV eff 
dr \ dr 



(68) 
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where 



y e// (0) = y(0) + -mV 



In view of the effective potential in (69), if one starts with 



one gets to 



V(4>) = -0 6 - 



(69) 



(70) 



(71) 



This is the <fi 6 potential - see Eq. (27); thus, it has the property that, if A is integer, the value m c = A shows that there 
will always be a critical m c for which there are analytical solutions representing topological solitons in the model. 
Such solutions are given by 



4>±{r)=± 



A 



„±2A 



(72) 



Moreover, for m < m c one can find non-topological solutions numerically; in Fig. [8] we show the potential for A = 2 
and for m = 0, 1, 2. The inset shows how the potentials behave for small: there are no lumplikc solutions for m = 
or 2, but for m = 1. 




FIG. 8: Plots of the potential in Eq. (71) for A — 2, and for m = 2,1,0, depicted as solid, dashed and dot-dashed lines, 
respectively. The inset details the behavior for (f> small. 



We now return to models described by a single real scalar field, and we consider the case D > 3. The equation of 
motion for <j) = <fi(r) is 

We write dx = ±r 1_D c?r to get d 2 (j>/dx 2 = WqW^q and again, we map the model into a one-dimensional problem. 
We solve dx = ±r 1 ~ D dr to get x — ipr( 2 ~ D ) /{D — 2). This shows that x is now in (— oo, 0] or [0, oo), and so we have 
to use the upper sign for x < 0, or the lower sign for x > 0. We can use the model of Eq. (38) for p — 4,6, ... to solve 
the D-dimensional problem with N = 2D — 2. In this case the solutions are 

^(Ap)( r ) =tan h p 



D-2 



(74) 
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for D = 3,4, ... Their energies are given by = £l£i[2p/(4p 2 — 1)]. 

In the case of D = 3 and p = 4 the zero mode is given by r]j^' 4 \r) = c^ 3,4 ' tanh 3 (l/4r)sech 2 (l/4r); we could not find 
the normalization factor explicitly in this case. The zero mode binds at the skin of the defect solution, concentrating 
around the radius R of the defect, which is given by R — l/4arctanh[(l/2) 1 / 4 ]. 

We further illustrate the case D = 2 with another very interesting model, described by the potential 



The equation of motion can be written as 



(75) 



r— r— = sin(0) cos(0) 
dr \ dr 1 



(76) 



or, using r(d/dr) = d/dy, 



The first-order equations are 



— T = sm(cj>) cos(0) 

ay 



dy 



± sin(^>) 



(77) 



(78) 



This is the sine-Gordon model, and the soliton solution has the form cf>(y) = ±2 arctan(e ±y ). If we turn back to the 
radial coordinate, we get 4>(r) — ±2 arctan(r ±1 ), which is depicted in Fig. [9]. The model maps an old solution, shown 
in Ref. [50]. It reproduces topological solitons in magnetic systems, as one can see in Ref. [51]. The same solution 
appears in the more recent work [5 2], in the form of delocalized vortex structure in the spin-flop phase in crystals with 
the C nv symmetry. Our solution shown in Fig. [9] maps the angle formed between the staggered vector I and the 
tetragonal axis, as it is shown in Fig. [12] of Ref. [52]. 




1.6 



0.8 



I r 

°' 2 4 (3 8^ 10 

FIG. 9: Plot of the radial solution cj>(r) = 2arctan(r), which represents a topological soliton for the model described by the 
potential of Eq. (75) . 

The solutions that we have found for D > 2 have a central core, and a skin which depends on the parameters that 
specify the potential of the model. They are stable, well distinct of other known global defects such as the bubbles 
formed from unstable domain walls[35-38]. Also, they are neutral structures, and may contribute to curve space time, 
and to affect cosmic evolution. They may become charged if charged bosons and/or fermions bind to them. For 
instance, if one couples Dirac fermions with the Yukawa coupling Y(<j>) = r 1_£> W^, the fermionic zero modes are 
similar to the bosonic zero modes that we have just presented. 



14 



We notice that we have solved the equation of motion using r 1_£1 'dr = ±dx. Although this identification works very 
naturally for the first order equations, we used the equation of motion to present a more general investigation, which 
may be extended to tachyons in D spatial dimensions. For instance, from the lumplike solution for D — 1 and p = 2 
we can obtain another lumplike solution <j)) 2 ' 2 \r) = (r— l) 2 /(r+ l) 2 , which is valid for D = 2. Another issue concerns 
the identification of the topological behavior of the above defect structures. We do this introducing the (generalized 
current-like) tensor 

It obeys d^J^ 1 ^ 2 '"^ — for i = 1, 2, • • • , D, which means that the quantities 

P»lt2-"»D-1 _ j0il»2-"*D-l ^gQ^ 

constitute a family of D distinct conserved (generalized charge) densities. We introduce the scalar quantity 

Pi = Pnii-io-J 1 * 2 - 40 - 1 = {-1) D {D - l)!(#/dr) 2 (81) 

which generalizes the standard result, obtained with j M = e^ v d v <j) - see the reasoning leading to Eq. (64). Thus, we 
define the topological charge as 

Q° = J drp 2 D = (-1) D (D - l)\fl D AW (82) 

which exposes the topological behavior of the new global defect structures that we have just found. 
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